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ABSTRACT 

We investigate the potential of second generation measurements of redshifted 21 cm radiation from 
before and during the epoch of reionization (EOR) to reconstruct the matter density fluctuations 
along the line of sight. To do so we generalize the quadratic methods developed for the Cosmic 
Microwave Background (CMB) to 21 cm fluctuations. We show that the three dimensional signal 
can be decomposed into a finite number of line of sight Fourier modes that contribute to the lensing 
reconstruction. Our formalism properly takes account of correlations along the line of sight and uses 
all the information contained in quadratic combinations of the signal. In comparison with the CMB, 
21 cm fluctuations have the disadvantage of a relatively scale invariant unlensed power spectrum which 
suppresses the lensing effect. The smallness of the lensing effect is compensated by using information 
from a range of observed redshifts. We estimate the size of experiments that are needed to measure 
this effect. With a square kilometer of collecting area and a maximal baseline of 3 km one can achieve 
lensing reconstruction noise levels an order of magnitude below CMB quadratic estimator constraints 
at L = 1000, and map the deflection field out to less then a tenth of a degree (L > 2000) within 
a season of observations on one field. Statistical lensing power spectrum detections will be possible 
to sub-arcminute scales, even with the limited sky coverage that currently conceived experiments 
have. One should be able to improve constraints on cosmological parameters by using this method. 
With larger collecting areas or longer observing times, one could probe arcminute scales of the lensing 
potential and thus individual clusters. We address the effect that foregrounds might have on lensing 
reconstruction with 21 cm fluctuations. 

Subject headings: cosmology: theory - diffuse radiation - large-scale structure 



1. INTRODUCTION 

The possibility of measuring redshifted 21 cm radiation 
originating before and during reionization has attracted 
a lot or interest recently, as it promises new insights into 
this poorly understood epoch, when the flrst radiation 
sources turned on and eventually ionized the universe. 
The origin of this observable is the spin flip transition of 
neutral hydrogen. Fluctuations in its brightness temper- 
ature can be caused by baryonic density fluctuations, in- 
homogeneities in the spin temperature, and during reion- 
ization also by an inhomogeneous ionization fraction. Fu- 
ture observations potentially contain information about 
all of these fluctuation sources (Barkana & Loeb 2005; 
Loeb & Zaldarriaga 2004). 

These observations could be the deflnitive probe of the 
'dark ages'. At present, there are several probes of the 
epoch of reionization. First there is the total Thomson 
scattering optical depth inferred from the CMB polariza- 
tion on large scales points to a mean redshift of reioniza- 
tion, Zri = 10.9^2 3 (Spergel et al. 2006). This constraint 
is prone to large additional systematic uncertainties, as 
we are just beginning to understand the properties of po- 
larized foreground emission (Page et al. 2006). The CMB 
offers only an integral constraint so in principle it allows 
for very complicated evolution scenarios of the ionized 
fraction, and does not contain much information as to 
whether the ionization process is homogeneous or not. 
More detailed polarization information will allow a few 
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more numbers characterizing the reionization process to 
be measured (Zaldarriaga 1997; Hu & Holder 2003). 

Quasar absorption spectra in the redshift range 5.5- 
6.5 indicate a fast evolution of the ionized fraction while 
xh > 10~^ (see e.g. Fan et al. 2005). They carry a Gunn- 
Peterson trough (Gunn & Peterson 1965), suggesting in- 
complete ionization of the intergalactic medium (IGM) 
along the sight line. The sizes of the HII regions directly 
surrounding the quasars can be used to infer a tighter 
constraint on the minimal neutral fraction, xh > 0.1 
since the IGM has to be somewhat neutral otherwise the 
regions would be larger (Wyithe et al. 2005). 

The redshifted 21 cm line is outstanding in comparison 
to CMB and quasar observation studies because it should 
be sensitive to order unity changes in the neutral fraction 
Xh, hence probing the middle stages of reionization. It 
will also provide well-localized measurements along the 
line of sight, and it will not require the background of 
bright sources which might be rare at high redshifts. 

Although 21 cm measurements offer the possibility to 
constrain the state of the IGM during the end phases 
of reionization, due to the confounding effects of their 
unknown astrophysics it might be difflcult to compete 
with the constraints on cosmological parameters (Mc- 
Quinn et al. 2005) coming from microwave background 
experiments such as the Planck satellite ^, the E and B 
experiment (EBEX) (Oxley et al. 2004). the Atacama 
Cosmology Telescope (ACT) ^ (Kosowsky 2003), or the 
South Pole Telescope (SPT) ^ (Ruhl et al. 2004)). 

* http: / /www.rssd.esa.int/index.php?project=PLANCK&page=perf_top 
^ see http://www.hep.upenn.edu/ angelica/act/act. html 
^ see http://astro.uchicago.edu/spt/ 



2 



In this paper we will explore the information contained 
in the intrinsieally three dimensional 21 cm measurement 
about the intervening mass distribution at lower redshifts 
through the lensing effect. We generalize the quadratic 
estimator technique developed by Hu (2001b) to a three 
dimensional observable. Our formalism uses all the in- 
formation in both shear and convergence and properly 
takes account of correlations along the line of sight. Be- 
fore the ionized fraction of the IGM becomes substantial, 
the 21 cm emission against the CMB is a near Gaussian 
random field. In this regime, which is the focus of the 
present article, the quadratic estimator should be close 
to optimal Hirata & Seljak (2002). 

We will look at 21 cm fluctuations at their lowest pos- 
sible redshift range, z = 6 — 12, where the second gen- 
eration of experiments might be able to measure at high 
signal-to-noise. Using our formalism we can also estimate 
information losses due to foreground contamination, once 
these are described by some model. We also explore the 
possibility to constrain the dark energy density with lens- 
ing of the 21 cm background. Another application would 
be to measure nonlinearities in the density field. We 
compare our results to the potential of a future high 
precision observation of the CMB. The CMB damping 
tail is an advantage for lensing reconstruction since re- 
construction errors decrease with increasing slope, but 
at the same time it leads to a small scale limitation for 
CMB reconstruction, as the signal quickly falls below the 
noise. 

Lensing reconstruction using the redshifted 21 cm ra- 
diation in absorption against the CMB has been inves- 
tigated by Cooray (2004), in particular also the possi- 
bility to get a handle on gravity waves from inflation 
by gravitational lensing cleaning of B mode polarization 
(lensing converts E to B modes of the polarization and 
acts as a contaminant to the primordial signal) (Sigurd- 
son & Cooray 2005). Their work describes two types of 
observations, which employ a 20/200 times larger total 
collecting area than we will assume here (and five times 
longer observation time) to observe angular fluctuations 
Lmax = 5000—10^ in the 21 cm field. The problem is that 
the prospect, beating the CMB level for B mode lensing 
cleaning, relies on measuring the 21 cm power spectrum 
at very high redshifts (they use Zgourcc = 30), at which 
the galactic synchrotron contamination is a factor ~ 20 
larger than for example at redshift 8. The reason such 
high redshift observations are needed to compete with 
likelihood based lensing estimation is that there is a par- 
tial delensing bias when comparing lenses out to different 
redshifts and Zg = 30 turns out to be close enough to the 
last scattering surface of the CMB, where gravity waves 
are expected to create the B mode fluctuations. These 
authors furthermore make the approximation of treat- 
ing their slices through the 21 cm measurement cube as 
uncorrelated, which is not warranted. 

Pen (2004) suggested measuring the effective conver- 
gence from the effect it has on the real space variance 
map of 21 cm fluctuations. The author also presented 
rough sensitivity estimates for LOFAR, PAST, and SKA. 
We improve on this by using all the available informa- 
tion in convergence and shear in an optimal way and use 
more realistic errors. 

In section 2 we will review the theory of the redshifted 
21 cm signal in different regimes, which leads to the dis- 



tinction of a neutral phase where the brightness temper- 
ature fully traces the baryonic fluctuations, as well as a 
patchy phase, which can be modeled separately. 

In section 3 we review the quadratic estimator tech- 
nique for lensing reconstruction following Hu (2001a). 
Then we naturally expand the formalism to the extrac- 
tion of weak lensing from an intrinsically three dimen- 
sional signal. The decomposition of the line of sight 
component of the signal into modes leads to a hierarchy 
of independent lensing backgrounds that can be probed 
with varying precision. Although our concentration lies 
on applying the quadratic estimator to the epoch imme- 
diately before substantial ionization occurs, we scrutinize 
its applicability to the patchy regime by using an analj^tic 
model for the morphology of HII regions. 

In section 4 we put our investigation in the context of 
experiments, estimating the potential redshift range in 
which they will observe. Because signal and noise are 
evolving with redshift, we break down the volume of the 
observation in smaller boxes along the line of sight. We 
calculate estimates for the lensing reconstruction future 
21 cm experiments might be able to achieve. 

We present results based on current rough specifica- 
tions of the Square Kilometer Array (SKA) in Section 5 
and compare them to the possibility of constraining the 
matter power spectrum with the CMB temperature and 
polarization. The 21 cm approach turns out to have no 
angular scale limitation for constraining the convergence. 
We address limitations due to the galactic foreground 
and show how these can be incorporated into our for- 
malism in a straightforward manner. We conclude with 
an outlook and discussion in section 6. 

A ACDM cosmology is assumed throughout all calcu- 
lations, with parameters il„i = 0.3, Oa = 0.7, Qt = 
0.04, Ho = lOOhkm/sec/Mpc (with h=0.7), and a scale- 
invariant primordial power spectrum with n = 1 normal- 
ized to (Tg = 0.9 at the present day. 

2. THEORY OF THE REDSHIFTED 21 CM SIGNAL 

The shape of the 21 cm fiuctuations is easy to predict 
in the absence of radiative sources or in a rather homo- 
geneous radiation background. Then it is simply pro- 
portional to the matter power spectrum which we know 
based on measurements of the CMB, galaxy surveys, etc. 
In the presence of fluctuations in the ionization fraction, 
the observable is more difficult to quantify, due to our 
lack of knowledge about the spectrum and efficiency of 
first sources, and the difficulty of modelling feedback pro- 
cesses. Some progress has been made using radiative 
transfer codes (Sokasian et al. 2003) on up to 20 Mpc/h 
scales, but because of the large size of the HII regions 
that simple analytic considerations suggest (Wyithe & 
Loeb 2004; Purlanetto et al. 2004a) it may be helpful for 
understanding the large scale morphology of HII bub- 
bles to also use analytic models, such as the one based 
on extended Press-Schechter theory by Purlanetto et al. 
(2004b). 

The optical depth of a region of the IGM in the hyper- 
fine transition is (Field 1959) 
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Here = 1420.4 MHz is the rest-frame hyperfine tran- 
sition frequency, Aio = 2.85 x 10~^^s~^ is the sponta- 
neous emission coefficient for the transition, Ts is the 
spin temperature of the IGM (i.e., the excitation temper- 
ature of the hyperfine transition), Tcmb = 2.73(1 -I- z)K 
is the CMB temperature at rcdshift z, and nui is the 
local neutral hydrogen density. In the second equality, 
we have assumed sufficiently high redshifts such that 

H{z) « Hoflll'^^l + z'f'l'^ . which is well-satisfied in the 
era we study, z > 6. The local baryon overdensity is 
\ ^ bs = p/p and xh is the neutral fraction. Here the 
index in the overdensity emphasizes that we are mea- 
suring quantities in redshift space. The radiative trans- 
fer equation in the Raylcigh- Jeans limit tells us that the 
brightness temperature of a patch of the sky (in its rest 
frame) is Tb = TcMBe"^ + Ts{l - e""). Then the ob- 
served brightness temperature increment between this 
patch, at an observed frequency u corresponding to a 
redshift 1 + z = vq/u, and the CMB is 

Ts — Tcmb 
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So if the excitation temperature in a region differs 
from that of the CMB, the region will appear in emission 
(if Ts > Tcmb) or absorption (if Tg < Tcmb) against the 
CMB. 

Ciardi & Madau (2003) argue that the intensity of 
the Lya background is large enough at all relevant red- 
shifts, so that Lya pumping (the Wouthuyscn-Ficld ef- 
fect (Wouthuysen 1952; Field 1958), in combination with 
the fact that the color temperature is in thermodynamic 
equilibrium with the gas kinetic temperature in an opti- 
cally thick medium, effectively couples the spin temper- 
ature to the kinetic temperature of the gas. Between the 
CMB last scattering at redshift z « 1000 and z « 140, 
Thomson scattering of the CMB photons effectively leads 
to Tcmb — Tk, but the gas decouples at 2 ~ 140 and 
cools adiabatically until significant structure begins to 
form. For the regime we are interested in, z < 15, it 
is likely that the first supernovae and/or accreting black 
holes produced enough high energy photons to heat the 
IGM Madau et al. (1996). In this case it can be assumed 
that a relatively homogeneous X-ray background leads to 
Ts « Tif > Tcmb (e.g. Ciardi & Madau 2003; Madau 
et al. 1996) so that the temperature factor in Equation 2 
becomes simply unity. In this picture inhomogeneous ef- 
fects on the brightness temperature due to shock heating 
along sheets and filaments are a minor contribution. 

The 21 cm power spectrum is measured in redshift 
space, hence it is is related to the actual density con- 
trast power spectrum through 

P,(k) = P,(fc)[l + //4]2 (3) 



where f is the linear growth rate, / = Tf-^j^- The cosine 
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of the angle between the line of sight z and the wavevec- 
tor k is denoted /Ufe. 



Here we will assume that in the regime of interest 
/ ~ 1, hence we obtain with T, ^ Tcmb that the power 
spectrum of the 21 cm brightness temperature fluctua- 
tion in the neutral regime is given by 
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In this paper, we shall be mainly concerned with 
quadratic estimator lensing reconstruction during this 
highly neutral regime. The extension to a patchy epoch 
is complicated by the presence of a connected four-point 
function contribution to the source field. On the level of 
the power spectrum this contribution acts as a sample 
variance term, correlating different k-modes. In the fol- 
lowing paragraph we shall discuss a semi-analytic model 
that can be used to describe the patchy regime. In the 
analysis section wc will briefly estimate the application 
of our estimator to this regime. 

The fluctuating 21 cm brightness temperature simply 
traces the matter fluctuations until the time when the 
IGM starts getting significantly ionized. Current radia- 
tive transfer codes do not allow dynamical ranges large 
enough to simulate the reionization morphology on scales 
above 20Mpc/h. To model the patchy regime we imple- 
ment the analytic model suggested by Furlanetto et al. 
(2004b) into the three dimensional signal Gaussian ran- 
dom fields of sidelength lOOOMpc/h comoving (compare 
Zahn et al. 2005) with a resolution of 512"^. In the ex- 
tended Press-Schechter formalism the condition for a re- 
gion to be ionized, C/coU. > 1 leads to a condition on the 
mean overdensity within a region of mass m (Furlanetto 
et al. 2004b), 

Sm > 4(m, z) = 5c{z) - ^/2K{0[ai^ - a\m)f/^ (5) 

where K{Q = erf~^(l — C~^); C is the ionization effi- 
ciency factor, a^{m) is the variance of density fluctu- 
ations on the scale m, and 6c{z) = 1.686/T'(a) is the 
critical density for collapse scaled to today, with D{z) 
the scale factor, a^^^ is the variance on the scale of fluc- 
tuations on the scale mmin the mass for which T = W^K 
at virialization. 

Figure 1 shows the angular power spectrum of 21 cm 
fluctuations at redshift 8 for different values of the effi- 
ciency and therefore ionization fraction. We see that up 
to an ionization fraction of Xi = 0.7 the signal on small 
scales decreases (structures below the bubble scale are 
washed out) while the ionized regions lead to a bump 
in the power spectrum on scales k ~ 0.2Mpc/h~^ which 
corresponds to an angular multipole at the relevant red- 
shifts of approximately I = 1000. At even higher ionized 
fractions, the entire patchy power spectrum lies below 
the neutral case, so it becomes more difficult to observe 
the fluctuations and hence to use them for lensing recon- 
struction. In any event it seems from Figure 1 that a 
significant part of the patchy regime can be used for the 
lensing reconstruction, in part because the boosted am- 
plitude on scales Z ~ 100 — 1000 might aid somewhat in 
this effort. We will evaluate these estimates more carc;- 
fuUy in Section 5. Our choice of ionization efficiency 
lets reionization begin at redshifts around 8 and evolves 
rather slowly, due to photon consumption by clumps in 
the IGM. It evolves through Xi ~ 50% at redshift 7 to 
complete at redshift 6. 




Fig. 1. — Power spectrum of the 21 cm signal at various stages 
during reionization. The thick solid line shows the expected noise 
for SKA for a flat configuration of antennas inside a circle of diam- 
eter 2 km. When the ionization fraction rises during the expansion 
of HII spheres, first the overall signal decays, then bubbles cluster 
quickly to form large ionized regions tens of Mpc/h across. This 
can lead to a significant increase in the signal on scales that the first 
imaging 21 cm experiments should be sensitive to. On the smallest 
scales the signal is decreased, making usage of 21 cm fluctuations 
for lensing reconstruction somewhat more difficult. 



We will assume a bandwidth ofB = 5MHz for calculat- 
ing the power spectrum at various redshifts. This corre- 
sponds to a redshift interval Az = 0.286. During the neu- 
tral phase the density fluctuations evolve slowly, however 
the sensitivity of the experiments change rapidly with ob- 
servation frequency. Because the comoving length scale 
is given in terms of bandwidth through 
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our window in frequency space corresponds to a depth 
of 50 Mpc/h at z=6 and to 70 Mpc/h at z=12. When 
the first extended HII regions start forming, the power 
spectrum evolves more rapidly, however reionization still 
only lapses over a comoving length of about 300 Mpc/h, 
in comparison to which our window is small. 

We will generate 21 cm brightness temperature power 
spectra following Formula 4 using linear power spectra 
that contain the acoustic oscillation amplitude generated 
from CMBFAST transfer functions (Seljak & Zaldarriaga 
1996). The baryonic wiggles included in the code might 
aid somewhat in our reconstruction endeavour, for rea- 
sons given in the next section. We use the same transfer 
functions to implement the signal in three dimensional 
Gaussian random fields and model the patchy phase. 

3. WEAK LENSING RECONSTRUCTION 

3.1. Quadratic Estimator, General Consideration 

Lensing by large scale structure happens whenever 
there is a fluctuating background field. At position n 
we observe the field 



where T denotes the unlensed field, 56 = \7(p is the dis- 
placement vector while (/) is the projected potential. Here 
and in what follows, boldface quantities denote vectors. 
The projected potential is given in terms of the gravi- 
tational potential ■ijj{x,D), where x is position and D is 
used as a time variable, as 
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So the power spectrum of the displacements is Cf^^^ = 
L{L + l)Ct^. 

An estimator D{n) for the lensing displacement field 

information contained in the temperature field should 
contain an even number of temperature terms since the 
expectation value for odd powers would vanish. It must 
also satisfy the condition 

{D{n)) = 60{n) (9) 

when averaged over many realizations of the background 
radiation field, that is for example the CMB or 21 cm 
radiation. 

Hu (2001a) showed that the divergence of the 
temperature- weighted gradient of the map achieves max- 
imal signal-to-noise among quadratic statistics. In 
Fourier space this quadratic estimator takes the form (Hu 
& Okamoto 2002) 

i)est.(L) = A(L) j ^i.(l,L-l)r(l)r(L-l) (10) 

The Filter F(l, L — 1) is obtained by minimizing the vari- 
ance of Dest. under the normalization condition for -D(L) 
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Here 0*°* is the sum of the lensed angular power spec- 
trum of 21 cm fiuctuations C; and the noise power spec- 
trum which we will give in the next section, Cj^ . As 
shown by Mandel & Zaldarriaga (2005), the effect of 
lensing on the angular 21 cm power spectrum for an 
individual plane is small, so one can use the unlensed 
power spectrum in place of it. The high-pass shape of the 
quadratic estimator gives it a property all reconstruction 
methods share, that they extract most information from 
the smallest scales resolved by some experiment. We em- 
phasize that the estimator is unbiased by construction 
(Equation 9), independent of Gaussianity of the source 
field. 

With the definition of the lensing reconstruction noise 
power spectrum Nl 

{D*{\.)D{V)) = (27r)2<5(L - L')(Cf ^ + N^) 

evaluation of the variance of Equation 10, ([[^(L) 
(277)2^^(0) ATf gives 
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Fig. 2. — CMB temperature and E polarization power spectrum 
and noise levels of our example experiments: Planck and a futur- 
istic polarized satellite mission (see text for specifications). The 
21 cm power spectrum has been rescaled to fit into the same plot. 
The CMB exhibits more pronounced peaks and a sharp decay on 
scales I ~ 2500 above which the Sunyaev-Zel'dovich effects become 
important. 



Roughly, when Li^Cf^^ = Nl- then structures down to 
the angular size Itx jL in the lensing field can be recon- 
structed. 

We illustrate the difi'erent power spectrum character- 
istics of CMB and 21 cm in Figure 2. 

We can investigate qualitatively which properties of 
the power spectrum set the effectiveness of this estimator 
by looking at the; limit L — > (on scales where the noise 
is small and the lensed power is not systematically larger 
than the unlensed, C^°* = C; + Cf ~ C; ~ C;) 
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The integrand of 



dVnl 

the first order term has its minimum close to a value a = 
—2, in other words when /^C; has no slope. This implies 
that the CMB with its exponential decay on small scales 
will have a lower value of Nl if those scales are resolved. 
Integrating Equation 16 to the maximum multipole Zmax 
at which Cf > C/^ leads to 
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where /min is the lower bound of the integration. The sec- 
ond term will lead to departures from a constant value of 
{L'^Nl). The scale at which both terms are comparable 
is 
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Evaluation of this expression shows that for a sloped 
power spectrum PCi, such as that of the CMB, L, 



comp. 



IS 



lower than for a spectrum that is nearly constant, such 

as that of 21 cm fluctuations. The second order term 
contributes strongly when the slope is large. If the mea- 
sured power spectrum has a small slope, as is the case 
with 21 cm fiuctuations, L'^Nl can be expected to be 
nearly constant to the scale Zmax where noise becomes 
important. 

In the following section we will generalize the quadratic 
estimator to a three dimensional observable that is used 
to reconstruct the lensing field. The final estimator is 
then simply the sum of estimators of k modes along the 
line of sight. This allows one to improve the constraints 
above those by the CMB especially on small scales by 
summing over sufficiently many lensing backgrounds. 

For our analysis in Section 4 we will compute the angu- 
lar power spectrum of the deflection angle as the integral 
over line of sight 
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over the power spectrum of mass fluctuations. Dj. denote 
angular diameter distances with x — l,s for lens and 
source. 

We use the halo model fitting function to numerical 
simulations of Smith et al. (2003) to generate nonlinear 
A CDM power spectra as input for the deflection angle 

integral. 

3.2. Extension to a three dimensional signal 

Different from the CMB, in the case of 21 cm bright- 
ness fluctuations we will be able to use multiple redshift 
information to constrain the intervening matter power 
spectrum. One could imagine applying this estimator to 
succesive planes prcpendicular to the line of sight. How- 
ever the different planes would be correlated and there is 
no straightforward criterion to take this correlation into 
account in establishing the final estimator, since it would 
depend on redshift, source properties, and during reion- 
ization on the mean bubble size and distribution. In- 
stead we will use the knowledge of the three-dimensional 
information in a different way, dividing the temperature 
fluctuations in Fourier space in fluctuations kj_ perpen- 
dicular to the line of sight, and a component fcy in the 
frequency direction. 

We devide a volume on the sky to be probed by a 
given 21 cm survey into a solid angle and a radial 
coordinate z. Components of wavevectors along the line 
of sight are described by k\\ and those perpendicular to 
the line of sight are given by the vector kj^, which is re- 
lated to the multipole numbers of the spherical harmonic 
decomposition on the sky as 1 = kj_r>, where T) is the 
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angular diameter distance to the volume element wc arc 
probing. 

Suppose we want to measure a field /(r) with power 
spectrum P(k) = P(kj_,fc||). Converting k± to angular 
multipole 
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with /(I, fc||) = we have 
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where in the second step a factor of V'^ got absorbed into 
^^(1 - 1') because of 1 = k_LD. 

Let us discretize the z direction of a real space observed 
volume with radial length £, 
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so that 
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SO that on the sphere we have 
(4(li)7;(l2)) = (27r2)5^(li-l2)5,,,, ] (27) 



where /x^ is the cosine between the wave vector and the 
line of sight. 

We have for the angular power spectrum for separate 
values of j that (including redshift space distortions) 



P2£ 



(28) 



where P now represents the spherically averaged power 
spectrum. We have also introduced the notation Cij, 
denoting the power in a mode with angular component 1 
and radial component kj = j 2t:/C. 

Wc show in the Appendix that because modes with 
different j can be considered independent, the best esti- 
mator can be obtained by combining the individual es- 
timators for separate j's without mixing them (in the 
sense of making quadratic combinations of them). As 
long as the IGM is not substantially ionized, the assump- 
tion of Gaussianity is justified at the redshifts of interest, 
2; ~ 6—12, where non linearites in the gravitational clus- 
tering are small on observed scales of several Mpc/h. We 
find the three dimensional lensing reconstruction noise 
defined by 

{D(L)D* (L')) = {2Trfd(L - L')iVf (29) 



to be (Equation A21 of the Appendix, where e.g. C 

1 



N 



D 



1 



(30) 



(31) 



where in the last line we have just substituted the stan- 
dard expression, Equation 15. Note that analogous to 
the CMB case, Cf* = C; + Cj^ , where Cf is the noise 
power spectrum, and C; ~ Ci in the case of 21 cm fluc- 
tuations. 

If there is a connected four point function contribution 
during the epoch of extended HII regions, this will add a 

term to the variance of the estimator. This will change 
the weightings Al as well, and make the estimator sub- 
optimal. We will not try to c;omc up with an estimator 
that is fully compatible with the non Gaussian signal due 
to patchy reionization. Instead, in the next section we 
will conservatively estimate that sensitivities get worse 
by a factor ^^,1— , where xh is the ionized fraction. In 

other words, we will effectively be treating the patchy 
regions as part of the source field which can be masked. 

Let us examine the contribution of the components 
Cij to the total noise in the estimation of the deflection 
field. From Equation 28 we see that with higher values 
of the component fc||, higher values of the three dimen- 
sional power spectrum P{k) will translate into each C; 
value. But P{k) is monotonically falling on all scales 
of interest (foregrounds can be expected to contaminate 
constraints below k < 10~^h/Mpc, the scale of the hori- 
zon at matter radiation equality and the turnover of the 
power spectrum), so that effectively the angular power 
spectrum amplitude will drop with j towards the noise 
level. We show this in Figure 3 for our bandwidth choice 
of 5MHz. Notice that on the smallest resolved scales, 
I ~ 2000 — 5000 the signal increases slightly as we go 
from the fundamental fc|| = to the next higher modes 
in our line-of-siglit decomposition. This is because the ra- 
dial component of the power spectrum is increased due to 
redshift space distortions. Because the smallest resolved 
angular scales contribute most to the lensing reconstruc- 
tion, in case of the first few modes in the fcy decomposi- 
tion this increase overcompensates for the general decay 
of Cij on the smallest scales. However overall for each 
rectangular data field of a given bandwidth we will only 
be able to sample a limited number of modes fcy with 
signal-to-noise greater than one, this number being pro- 
portional to the frequency depth/bandwidth of the field. 

The 21 cm power spectrum and, as we will sec in the 
next section, the noise of the experiment depend on red- 
shift, so we calculate them for volume elements corre- 
sponding to the above choice a 5MHz bandwidth. Be- 
ginning at the end of reionization, each volume element 
contributes to reducing the reconstruction noise while 
wc go along, until the signal-to-noisc of the experiment 
at high redshift becomes negligible. \\v should wonder 
whether long wavelength modes overlapping neighbour- 
ing volumes lead to an underestimation of the final lens- 
ing noise, as we are treating each volume element as un- 
correlated. We made a simple test by comparing the Nl 
of a pair of neighbouring volumes to the sum of A''i:,'s 
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Fig. 3. — The dotted lines show the angular power spectrum for 
different fcy ^ = j2iT/L, labeled by j. Going to higher values of 
jr, the signal decays quickly below S/N=l and only the first 15 or 
so modes contribute to the final lensing estimator. This number 
depends of course on the thickness of the redshift interval probed 
(in this case B = 5MHz, corresponding to Az = 0.286 at z = 8). 
T?he thick and thin solid curves are for a flat and a cored antenna 
configuration of SKA (see the text). 



of each clement. With our choisc of a 5MHz bandwidth 
we found no excess information > 1% in the stim of the 
separated field's A'/^'s, meaning that wc; can safely neglect 
those correlations. A multiplicity of those volumes can 
thus be used to reconstruct each lens just by summing 
over them. 

4. ANTENNA CONFIGURATION, SENSITIVITY 
CALCULATION 

We have found out in the previous section that in or- 
der to measure the lensing signal with redshifted 21 cm 
fluctuations, we need small angular scale resolution as 
well as wide redshift coverage. Sensitivity calculations 
for 21 cm experiments are presented in Morales (2005) 
and McQuinn et al. (2005). The sensitivity diminishes 
quickly with rising foreground temperature at longer 
wavelengths. The foreground temperature is dominated 
by galactic synchrotron. To observe at high redshift one 
needs to compensate for the increased foreground by in- 
creasing the collecting area or number of antennas. The 
angular resolution is improved by increasing the maximal 
baseline. To have good angular Fourier mode coverage 
we want all baselines to be represented though. Hence 
an increase in angular resolution has to be compensated 
for by increasing the collecting area in order to keep the 
covering fraction the same. In this section we will dis- 
cuss various aspects of the first two generations of 21 
cm experiments in the context of using them for lensing 
reconstruction. 

A standard approach in radio astronomical measure- 
ments (see e.g. the Very Large Array (VLA) , the 
Low Frequency Array LOFAR ^ , and the Atacama large 
millimeter array (ALMA) ^ configurations) is to have 

http://www.vla.nrao.edu/ 
* http;//www.lofar.org/ 
® http://www.alma.nrax).edu/ 



a power law decay in the number density of antenna 
V cc r" with radius. This translates into a drop in the 
number of baselines with separation. The distribution 
flattens out towards the center to n{r) = r° simply be- 
cause antennas cannot be stacked closer to each other 
than their individual physical size. We found in the pre- 
vious section that the contribution to the weak lensing 
estimator of Cij becomes smaller quickly with higher line 
of sight modes j used. This is just an expression of the 
fact that what we are interested in is displacements of 
21 cm photons perpendicular the line of sight. Lensing 
reconstruction works through a large number of individ- 
ual sources being aligned around a big deflector, hence 
the imaging of small scales is important. For probing 
small angular scales a flat array profile (that is without 
a power law decay) may offer an advantage, especially if 
the observable angular power spectrum falls slightly be- 
yond I ~ 1000 and we want to probe it on those scales. 

We can calculate the number of baselines as a function 
of visibility u from performing the autocorrelation 



n{u) j d^l' Aground (r + X) 'Pground(r) 



(32) 



where A is the observed wavelength and V{r) denotes the 
radial profile of the circularly symetric antenna distribu- 
tion. X is the vector of seperation of an antenna pair, 
X = Au. 

The time a particular visibility u is observed, tu is given 

by 

tu = ^n{u) (33) 

where is the effective antenna area and to is the to- 
tal observing time. We will assume 2000 hours for our 
calculations, which might be achieveable with planned 
observatories within a single seasons. 

The sensitivity for a given array distribution and spec- 
ifications can then be calculated as follows. The RMS 
fluctuation of the thermal noise per pixel of an antenna 
pair is (White et al. 1999; Zaldarriaga et al. 2004) 



AeVBt 



(34) 



where B is the bandwidth of the observation. For a single 
baseline the thermal noise covariance matrix becomes 



a 



N 



"1-3 

Bt^ 



(35) 



where T^ys is the system temperature, dominated by 
galactic synchrotron radiation (roughly, T^ys oc I/-2-55 
with Tgys = 440K at z = 6) and B is the bandwidth of 
this frequency bin of the total observation. From this we 
get the noise versus angular multipole number because 
of 2iTU = I through |^C/^ = d'^uC^ . 

We are now going to asses the potential of planned 
21 cm experiments to measure the lensing imprint using 
this quadratic estimator. There are currently four ma- 
jor experiments under way, the Mileura Wide Field ar- 
ray (MWA) (Morales 2005), the primeaval structure tele- 
scope (Pen et al. 2004), the Low Frequency Array (LO- 
FAR) ^° and in the second generation planning stage, the 

http://www.lofar.org 
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Square Kilometer Array (SKA) While these projects 
differ qualitatively in the hardware used, the most cru- 
cial difference is their total collecting area. For SKA, 
current proposals call for 5000 antennas with an indi- 
vidual collecting area of 120m^ at 2; = 8 (the effective 
antenna area depends on wavelength and hence redshift 
through the square). The collecting areas for MWA and 
LOFAR/PAST are significantly smaller, 1% and 10% of 
that of SKA respectively. The optimal array can be de- 
signed by distributing the total collecting area over a 
large number of antennas. The number of baselines goes 
as -/V^„( which enters into Equation 33. In other words, 
since an array with say ten times the number of dishes 
(each ten times smaller) as SKA has a ten times larger 
survey speed, it is better suited for EOR observations. 
However the computational requirements for the associ- 
ated correlator unfortunately also scale as N^^j^i ■ 

We show the result of our imaging sensitivity calcula- 
tion for a flat (r°) to r^ax = 1500 m and a cored (r" to 
r = 80 m, then r^^ to Tmax = 1500 m) array of SKA an- 
tenna in Figure 3, together with the hierachy of decreas- 
ing angular power spectra as the line of sight component 
of the 3D measurement is decreased. We find that a 
constant radial density of antennas turns out to offer the 
best compromise between the number of (kj_, fc||) modes 
that can be probed and overall angular resolution. For a 
cored r^^ distribution for example, the improved sensi- 
tivity on those intermediate scales does not compensate 
for the loss in angular resolution that comes from the 
lack of large baselines relative to a flat distribution. 

For statistical detections of the convergence on the sky 
a large field of view is desired. The angular power spec- 
trum errors go as 1 / -y/Jsky • Hence it is advantageous to 
distribute the same total collecting area in many smaller 
antennas, thus increasing the speed of the survey. On the 
down-side, it is more difficult to scale small dipole like 
antenna to high freqencies (below the redshift of reioniza- 
tion), where other interesting science with 21 cm emitters 
lies (see e.g. Zhang & Pen 2005). 

5. RESULTS AND A COMPARISON WITH THE CMB 

In this section wc will establish a specific configuration 
that seems likely within reach of the second generations 
of EOR observatories and calculate our three dimensional 
lensing estimator for a particular reionization scenario. 
This way we will gauge the prospects of lensing recon- 
struction using redshifted 21 cm fluctuations and how 
they relate to using other backgrounds. 

The prospects of lensing reconstruction with the 
quadratic estimator have been explored in depth by Hu 
(2001a); Hu & Okamoto (2002); Hu (2002); Amblard 
et al. (2004) in the context of the CMB. A systematic 
problem with temperature is that it is contaminated 
by Dopplcr related anisotropics (the kinetic Sunyaev- 
Zel'dovich effect Sunyaev & Zeldovich 1980) which share 
the same frequency dependence as the primordial CMB. 
Amblard ct al. (2004) showed that this contamination is 
significant, and might eventually have to be taken care 
of by masking thermal SZ detected areas from the map. 
It is found that using the CMB polarization an order 
of magnitude improvement over temperature reconstruc- 
tion might be possible (Hu & Okamoto 2002). On scales 

http://www.skatelescope.org/ 



where the B component of polarization can be resolved 
(above I ~ 100 this component is almost entirely pro- 
duced by gravitational lensing rotation of the E modes), 
with iterative likelihood techniques astonishingly there 
may be no limit at all to delensing (Seljak & Hirata 2004), 
on angular scales where the B modes can be resolved. 

As discussed in the previous section (Equation 16), 
the effectiveness of lensing reconstruction depends on 
the slope of the power spectrum, and the 21 cm back- 
ground suffers a disadvantage compared to the CMB in 
that there is no damping tail and the traces of baryonic 
oscillations arc comparatively small. On the other hand 
the lack of an exponential decay in the 21 cm fluctua- 
tions suggests that these can be used for reconstruction 
out to smaller scales. Figure 2 shows CMB strengths 
and limitations in comparison to the shape of the 21 cm 
power spectrum (which has been rescaled to fit in the 
same plot). 

The damping tail of the CMB at scales I > 3000 leads 
to a limit for the deflection scales I ~ 1200 that can 
be probed using quadratic techniques (Hu & Okamoto 
2002). The 21 cm power spectrum docs not decrease 
substantially at small scales, the only theoretical limit 
being set by the Jeans scale. 

Figure 4 shows a rather constant shape of L'^Nl ver- 
sus the angular multipole of the lensing field L. This is 
essentially a consequence of the rather fiat shape of the 
21 cm power spectrum, as we discussed following Equa- 
tion 16: in this case we would expect Nl to not change 
much up to scales where the signal becomes comparable 
to the noise of the experiment. This property will allow 
us to image the lensing field down to small scales. Wc 
also show in the same figure a plot of the temperature 
based CMB lensing reconstruction based on the Planck 
satellite experimental specifications. 

The CMB power spectrum sensitivity is given in terms 
of detector noise and beam (Tfwhm (Knox 1996) 



^noise _ y,-lgK'+l)<^FWHM/8'n2 



(36) 



We assume two types of experiments, Planck, with 
w~^^'^ = 27/LtK — arcmin at Sarcmin angular resolution, 
and a futuristic experiment with noise level w^^^^ = 
3/iK — arcmin at an angular resolution of 3 arcmin. For 

1/2 / — 1/2 

the polarization we use as usual that Wp = ^2wrp 

if all detectors are polarized. We include the noise power 
spectra for Planck temperature, and for polarization and 
temperature measurements of our reference CMB exper- 
iment in Figure 2. 

We calculate the minimum variance lensing noise level 
foUowing Hu & Okamoto (2002) 

^-(L) - ^ (37) 

where a and f3 run over T, E, and B temperature and 
polarization fluctuations. The noise levels N{'L)afs take 
on different forms depending on whether a and are 
equal (the BB term is generally considered to have van- 
ishingly small signal to noise), or different, a(3 = 6E, 6B, 
and EB. 

We compare Planck constraints to 21 cm observation 

sensitivities for the deflection angle power spectrum that 
we get from observing a 5MHz slice at redshift 8 in Figure 
4. The corresponding redshift depth is Az = 0.5. If we 
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Fig. 4. — Lensing reconstruction noise Ni^ for one redshift inter- 
val centered at z = 8 corresponding to a bandwidth of 5MHz. The 
curves labeled 'z=8' and 'z=1089' are the displacement field power 
spectra for 21 cm and CMB as source respectively. The thick solid 
curves labeled 'Planck' and '21 cm' are the lensing reconstruction 
noises we find. The 21 cm noise is based on combining all fcy modes. 
The thin lines labeled 'j=0-20' on the other side are the results for 
individual A;|| -modes. We see that from this redshift range alone the 
combined temperature and polarization information of Planck can 
be beaten. The necessity to substract foregrounds lessens the con- 
straint from 21 cm somewhat. They effectively render the first few 
fey modes useless for the reconstruction, see the text. The resulting 
noise levels axe shown in the thick dashed curves labeled jmin = 3 
and :?'min = 10, for a less and more conservative assumption about 
the complexity of foreground contamination respectively. 



want to gather all the information accessible to us from 
fluctuations in the 21 cm brightness temperature, we can 

add up a number of redshift intervals with the constant 
bandwidth, employing the whole redshift range covered 
by the observation. We find that even when using only 
a fraction of the entire data volume available in 21 cm, 
A2: ~ 0.3, Planck's combined temperature and polar- 
ization potential for doing lensing reconstruction can be 
beaten with the three dimensional generalization of the 
quadratic estimator. If we want to use a larger redshift 
range for the 21 cm reconstruction, we need to take into 
account the redshift evolution in the power spectrum of 
matter fluctuations, and that of the ionization fraction. 
For an SKA type experiment in the optimized configu- 
ration we use, this range is z — 6 — 12, at higher red- 
shifts the sensitivity becomes too small for imaging, the 
limitation being set by the increased temperature of the 
foreground. Note that if reionization completes earlier 
than z = 12 lensing reconstruction using 21 cm will be 
impossible with SKA. In our particular reionization sce- 
nario we achieve a constraint at L ~ 1000 about ten 
times as good (i.e. an Nl ten times lower) for the lens- 
ing reconstruction, in comparison to our reference CMB 
experiment with 3 arcminute resolution and a noise level 
of w~^/'^ = 3/iK — arcmin. We show this result in Fig- 



ure 5, together with nonlinear/linear lensing field power 
spectra. The Figure also shows (in the thick dashed 
curve) the increase in noise when our patchy model is 
used for the range 6-8 in place of an extension of the neu- 
tral phase. This increased noise relative to the case of 
probing the pre-reionization IGM has two reasons: first, 
as shown in Figure 1, the fluctuation level of the 21cm 
signal is decreased on the smallest resolved scales, mak- 
ing this source less valuable for lensing reconstruction. 
Secondly, the connected four-point function contribution 
adds a non Gaussian term to the noise covariance ma- 
trix of the power spectrum, acting as a sample variance 
term in correlating different band-powers. We treated 
this term in a simplified manner by assuming that the 
majority of bubble features arises on scales well above 
the resolution scale of SKA, about 1'. Then the bubbles 
can be resolved and masked when establishing the final 
estimator. The mask increases the sample variance sim- 
ply as 5Cf = —^SCi, where xh is the neutral fraction 

at this redshift. Our estimate of the lensing reconstruc- 
tion noise for the patchy epoch is conservative, in that 
power in the regions outside the bubble mask will not 
be suppressed, however we use the global average power 
spectrum which is suppressed by xjj. For completeness 
we also show (in the thick dotted curve) the total lens- 
ing reconstruction noise if we only use the neutral regime 
above 2; = 8 in the analysis. 

Because of the multiple background information, in 
reconstructing the deflection fleld we can approach the 
limit imposed by the finite angular resolution of an exper- 
iment (determined by its maximum baseline). In Figure 
6 wc show contours for the maximum lensing defiection 
field multipole L that can be probed as a function of the 
maximum baseline and the total collecting area. For the 
r'^ array there exists an optimal rmax which depends on 
the amplitude of the power spectrum and the total col- 
lecting area. If one uses a wider radius while keeping 
the number of dishes constant, the noise level increases 
on all angular scales, so although one might be able to 
use slightly larger multipoles (smaller scales) for the re- 
construction, many modes of our hierarchy do not enter 
the estimator. On the other hand, if one decreases the 
array size, thus making the antenna distribution denser, 
this leads to a lower noise level on relatively large scales, 
while no baselines exist anymore to measure the small 
scales that are crucial for the reconstruction. 

We find that after gathering 2000 hours of data, a very 
ambitious experiment with four times the collecting area 
of SKA (for which ^coU.,z=8 = 0.6km^) could detect lens- 
ing at scales beyond L ~ 6000, which is comparable to 
the scale of galaxy clusters. 

We will extend considerations to statistical detections 
of the deflection power spectrum. Because the geometri- 
cal ratio is larger for the CMB (qualitatively the 
bulk of lensing happens at angular diameter distances 
that are closer to the middle between the observer and 
the CMB), the respective defiection angle power spec- 
trum has a higher amplitude. On the other hand the 
21 cm experiment will have the advantage of measuring 
multiple planes so that at comparable angular resolution 
smaller errors in the angular power spectrum at high L 
can in principle be achieved. 

In Sigurdson & Cooray (2005) it was proposed that 21 
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Fig. 7. — Plot of the parameter F defined in the text. This 
suggests that 21 cm fluctuations are better suited in principle to 
measure the value of the cosmological constant than the CMB, 
assuming that the same number of angular modes can be probed. 



Fig. 5. — Lensing reconstruction noise with an experiment that 
has the total collecting area of SKA, if the IGM is not ionized 
during the regime z=6-12, shown in the thick line at bottom. The 
individual contributions from redshift intervalls are shown in the 
dashed lines. The noise levels are compared to the deflection angle 
power spectrum, where the solid (dotted) lines are for the nonlinear 
(linear) density fluctuations at z=8 and z=1089. On scales of L = 
1000 and above, our method might be able to achieve an order 
of magnitude lower total noise levels than what is possible with 
the CMB quadratic estimator technique (shown in the other thick 
solid line). Here the smallest angular scale reconstructed using 
EOR fluctuations is Lmax = 2250. The thick dashed line shows 
the increased noise level if reionization is inhomogcncous during 
2 = 6 — 8, calculated in the way detailed in the text. Finally the 
thick dotted curve gives the total noise level if the regime z = 6 — 8 
is not used at all in the analysis. 




AcoU.,z=8 = Na„tAe (km) 



Fig. 6. — Contour plot of Lmax, the largest displacement field 
multipole probed. From this Figure one can infer the optimal array 
radius (assuming here an distribution of antennas) that one 
should choose given a total collecting area A^oii ■ It also shows that 
if one would have four times the collecting area of SKA (for which 
Acoii. z=8 = 0-6km^), or if one were to observe four times longer on 
the same patch of the sky, very high multipoles could be probed. 



cm reconstruction of the lensing power spectrum might 
be helpful in getting at B mode polarization from pri- 
mordial gravity waves. Since lensing partially converts E 
(gradient) polarization into a curl component, this sec- 
ondary signal swamps the B modes from a possible in- 
flationary tensorial fluctuation background. The prob- 
lem is that at significantly lower redshifts than the last 
scattering surface, only part of the lensing structure en- 
countered by the CMB photons is traced, leading to a 
delensing bias if 21 cm fluctuations are used. One would 
have to observe at high enough redshifts if one were to 
compete with CMB polarization experiments. Indeed the 
authors find that in principle with an ultra sensitive ex- 
periment (for instance space based) observing at redshift 
30 one might be able to beat CMB limits beyond the 
iterative likelihood approach of Seljak & Hirata (2004); 
Hirata & Seljak (2003). To achieve this one would need 
1-2 orders of magnitude more collecting area than what is 
planned for the second generation of observatories, hence 
this application is beyond the scope of our paper. 

A characteristic of 21 cm fluctuations is that the dis- 
tance between observer and lens is a larger fraction of 
the whole distance to the source, so one would imagine 
that if the same number of multipoles are probed with 21 
cm as with the CMB (by having a lower noise level), one 
would obtain different constraints on dark energy. We 
show this in Figure 7, where the ratio 



rQ.(i) = 



(38) 



is plotted against L. The quantity measured shows how 

wcill the dark energy density can be constrained by using 
information from a given angular scale. Its value gives 
the constraint that can be put on the parameter from this 
scale if there were no degeneracies with other unknowns. 
We see that when the same range of angular scales are 
resolved, the 21 cm fluctuations fair somewhat better in 
constraining the dark energy density f^A- 
The noise for estimation of bandpowers is reduced by 
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averaging over L directions in a band of width AL 

A comparison between the polarized reference CMB 
experiment and an experiment with SKA's sensitivity is 
shown in Figure 8 where we assumed a sky coverage of 0.8 
for the CMB experiment and a smaller field of view 0.08 
for the 21 cm experiment. We plot the sample variance 
error ('S') and sample variance plus noise ('S+N'). Polar- 
ized CMB experiments suffer from foreground contami- 
nation (Page et al. 2006), so our reference experiment 
should be an idealized limit. The sky coverage of the 
21cm experiment might be achievable within one year of 
observation with an MWA type experiment that has the 
same collecting area as SKA but a ten times higher sur- 
vey speed. Using fluctuations in the 21 cm background, 
we should be able to measure the deflection power spec- 
trum to much higher multipoles, I > 10000 for this type 
of experiment, than what will ultimately be possible us- 
ing the CMB quadratic estimator technique, even if the 
latter observes on a much larger part of the sky (notice 
also that 21 cm experiments should be able to observe 
separate fields in consecutive seasons). The requirement 
of a combination of large collecting area and high survey 
speed makes this an ambitious project. 

Finally we would like to estimate the effect foregrounds 
will have on the estimation of the lensing field suggested 
here. As pointed out by Madau et al. (1996); Zaldar- 
riaga et al. (2004), fluctuations in the gas at high redshift 
can be scperated against the much brighter foregrounds 
(the main source of confusion being galactic and extra- 
galactic synchrotron), because the former vary rapidly, 
the latter slowly in frequency. This can be done by fitting 
various smooth functions to the signal. Morales et al. 
(2005), Wang et al. (2005), and McQuinn et al. (2005) 
suggest quadratic or cubic polynomials, but also more 
complicated functions such as Chebychcw polynomials 
have been proposed. McQuinn et al. (2005) show that 
this will practically make the first few wavevector modes 
in the line of sight direction unusable. The exact num- 
ber of modes that can be used depends of course on the 
nature of foregrounds, the bandwidth, and the technique 
of fitting. 

On the positive side once a model for foregrounds is 

given, wc can implement this pretty straightforwardly 
within our formalism by discarding the first few modes 
A;|[ (depending on the order of fitting polynomial) in the 
hierarchy of Cj's, compare Figure 3. It is to be expected 
that a large number of modes still will contribute to the 
final noise level. This is quantified in Figure 4 where 
the sum N^^ = V Efc,>,_^._^ ^i,fej£/(277) is plotted for 
the two cases jmin = 3,10 meaning that the first 3/10 
modes have been discarded in the analysis. The result- 
ing noise levels (short dashed lines) are somewhat higher 
than the solid line that was obtained from assuming no 
foregrounds. 



6. DISCUSSION 

In this paper we have extended the quadratic estima- 
tor formalism to tisc> a thrcK^ dimensional signal as lensing 
background. 21 cm fluctuations from neutral hydrogen 
prior and during the epoch of reionization contain an 
enormous amount of data points. The correlations in- 
duced by lensing into this signal can be used to probe 
the intervening matter fluctuations either on a individual 
object basis, or statistically to for example probe dark 
energy models. To describe fluctuations in the neutral 
fraction, we used an analytic model for the morphology 
of HII regions to demonstrate the applicability of our 
method to this regime. 

Our estimator should be complete as long as non lin- 
earities in the signal (e.g. due to the bubbles) are small. 
The bulk of the information we use is coming from the 
neutral phase. 

The first generation of experiments is likely not going 
to be able to image the angular fluctuations in 21 cm 
needed to measure the lensing effect. However we arrive 
at good constraints by employing current specifications 
of the SKA with a fiat antenna distribution. 

In comparison with the CMB, 21 cm fiuctuations have 
a rather featureless scale invariant unlensed spectrum 
which leads to a smaller lensing effect. This can be 
compensated for by using multiple redshift information 
when probing each individual lens surface. The CMB 
quadratic estimator sensitivity can actually be beaten 
this way, for example with an SKA type experiment. 

Another possibility would be to combine 21 cm lensing 
reconstruction with that from other observables, such as 
the CMB. Similar to combining the latter with galaxy 
shear surveys, this will improve the constraints on the 
mass/energy budget or geometry of the universe signifi- 
cantly. 

If ambitions in the community of observational and 
theoretical cosmologists increase in the years to come, 
and third generation experiments will be scheduled, a 
new prospect would be the measurement of polarized 21 
cm emission. Babich & Locb (2005) find that Thomson 
scattering of the quadrupole produced by the reionized 
universe produces the largest effect and similar to the 
gains by using CMB polarization, this could make the 
method suggested in this paper, the three dimensional 
form of the quadratic estimator, even more promising. 

Larger collecting areas and longer observation times 
also promise to allow 21 cm lensing reconstruction to 
map out the gravitational potential of individual galaxy 
clusters. We will attempt to address this topic in a future 
work. 
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Fig. 8. — Displacement field power spectra and sample variance ('S') and noise ('N') errors. An SKA like experiment might be suited 
to probe a large dynamical range in the displacement field. The displacement field angular power spectrum with error bars on the left is 
for our reference CMB experiment at /gjjy = 0.8 and the curves on the right are for 21 cm redshift 6-12 lensing reconstruction using the 
"MWA50k" (same collecting area but 10 times higher survey speed than SKA) with /gj^y = 0.08. 



APPENDIX 

QUADRATIC ESTIMATOR APPLIED TO A THREE DIMENSIONAL OBSERVABLE 
We want to observe the three dimensional field 1(0, z). We assume weak lensing 

1(61, z) = 1(61, z) + 5e- Vel{0, z) + ... 

where 1{9, z) is the lensed, 1{9, z) the unlensed field. The Fourier transform of this expression is 

dH' 



(Al) 



1(1, A:) -I 



;{5e{\-\')-i\)l{\',k) 



(A2) 
(A3) 



(2^)2 

= 1(1, k)~ j d?l'l(V, fc)0(l - 1') (1 - 1') • 1' 
where we have used that 59{L) — iL(f){L). We are looking for a quadratic estimator <I>(L) for (/)(L), i.e. of the form 

^JlSrJ^J If ^^^^^^ " ^' 

(notice that F{1, ki, k2, L) = F{1, k2, fci, L)). Because (5$(L) = (5$*(— L) it can also be shown that 

F(l,fci,fe,L) -F*(-l,-fci,-fc2,-L) (A5) 

We want to find F such that it minimizes the variance of $(L) under the condition that its ensemble average recovers 
the lensing field, ($(L))i = 0(L). This becomes (to first order in 



X (L - 1 - 1') • 1') - (I(L - 1, k2) J ^ 1(1', fci)0(l - l')(l - 1') • 1') 

J (2^)2 J 2tt J 2n ^ ' 1' 2' 



-l-l')x 



(271)^5(1 + l')(27r)(5(fci+A:2)P,,fcX 



X0(L - 1 - l')(L - 1 _ 1') • r + y -0^ (2^)25(L - 1 + l')(27r),5(fci + fe)Pi_,,fc0(l - l')(l - 1') • 1' 
'IS^J^J ^^afci'^2,L)(27r)<5(fci+fc2)[P,fc<^(L)L.(-l)+P^_,,0(L)L.(-(L-l^^^ 



/ (S^ / ff / ffm,fci,^2,L)(2^)<5(fci+A;2)[P,,fc0(L)L.l + Pi_,,fe0(L)L.(L-l)] , 



(A6) 
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where e.g. P;,fc is the power in a mode with angular component I and radial component k. With the requirement that 
($(L))i = <^(L) this leads to the normalization condition 

/ J^J ^ + ^2) [P/,feL • 1 + PL-;,feL • (L - 1)] = 1 (A7) 

The next step is to minimize the variance 

/llI/,Ml9^ f d^l f dki f dko f cPl ' f dki' f dko',, , , ^ ,n 

(2^7 277 (2^0^ y^y^F(U„.„L)F*(l',fci,.^,L')x 

X (1(1, fci)i(L - 1, fc2)i(l', k[)i(L - 1', k'2)) 

= lljjj J F{l,kuk2,L)F^l\k[,k'^,V){2w)H{l-l'){2n)6{k,-k[)P'i%x 

x{2nfS{l' -l){2n)6{k2- k',)Pr_i^,^ 

+ 11111 1 Fihki,k2,L)F%V,k[,k'^,V){2n)H{l-L + l'){2n)d{k^-k'^)P'i%x 

= J (0)^1^1 $(2-)^'5(o)F(l,fcl,fc,,L)^^*(l^fc;,fc^,L0PttP"-^,.. 

+ |f(2.)^^(0)i^(Ul,fc2,L)^^*(L-l^fc^,fc;,LOP^,\Pi°l,,, (A8) 

but from A4 we sec with the substitution L — 1 ^ 1 that F*(L — 1, A:2, fci, L) = F*{1, ki,k2, L) hence 

mL)\f) = 2(27r)2<5(0) J^J^J ^F{l,k„k„L)F*{l,k„k„L)Pl%Pl% . (A9) 

Both real and imaginary part of | |F| p = F'^+Ff contribute to this variance, however the condition for the minimization 
will only pick out the real part. The solution is found by minimizing the function 

(||$(L)||2) -ArX (EquationA?) (AlO) 

with respect to F, whore A;? is a Lagrangian multiplier. In steps, 

and 

djlMi^W) _o^o ^2^.n^ d^l dki dk2 , x ^ptot ptot .AIOA 

5i^(l,fci,fe,L)-2^2^) 5(0)^ — — 2FH(l,fci,fc2,L)P,,,^P^_,,,, (A12) 

so 

Fh(1, fcr, L) = An{2n)6'^{ki + fc.) P''^^^ ' ' + f ' " (A13) 

and by inserting this into the normalization condition A7 we get that 

r d^i [P,,feL.l+Pz._,.fcL.(L-l)]^ • (^1^) 



By using Equation A13 we find that the variance becomes 

{mL)f)=2{2nf6{0) f^J^J ^Pf%Pi°l,,,F| (A15) 

[P,,feL-l + Pi_,,feL-(L-l)]^ 



Using (2jr)J"(0) = g and / g|f ^ El,, "lis becomes 



ptot ptot 
^i,fe '^L-Lk 



(A16) 



(||$(i)||2)=2(27r)2^(0)A|j5] 



d^? [P;,fcL-l + P^_;,fcL.(L-l)]^ 

(27r)2 ptirt ptot 



k ^ ' l,k L-l,k 

= 2{2nr6ml^ 
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With the definition of the noise power spectrum iV*, 

($(L)$* (L')) = (2^)'<5^(L - L')iVf (A18) 
or(||$(L)|H = (27r)2(5^(0)iV* (A19) 

it follows that ^ 

" >p r d^i [P,,,L l-rP^_,.,L (L-l)]^ (^20) 
l^k J (27r)2 2Pt°* P*.°* , ^. 

and since the variances of estimator D of the displacement (-D)ens = and $ are just related by (||£)(L)|p) = 
L2(||$(L)||2) we obtain 

D 1 



V- 1 f d^i [Pj.fcL-l+Pj. i,fcL-(L-l)]^ 
Z-(fe L^' J (2-7r)^ 2PJ°.'P*-°* 

(A21) 



1 



so we have justified that we can simply sum over seperate k modes to arrive at our final lensing noise. 
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